Interference of stochastic resonances; splitting of Kramers' rate 
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We consider the escape of particles located in the middle well of a symmetric triple well potential 
driven sinusoidally by two forces such that the potential wells roll as in stochastic resonance and the 
height of the potential barrier oscillates symmetrically about a mean as in resonant activation. It 
has been shown that depending on their phase difference the application of these two synchronized 
signals may lead to a splitting of time averaged Kramers' escape rate and a preferential product 
distribution in a parallel chemical reaction in the steady state. 
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The escape of a particle from a metastable state due 
to thermal activation has been a major issue in chemi- 
cal dynamics and condensed matter physics for several 
decades ^[5]. As a typical paradigm in this context, con- 
sider a Brownian particle in the middle well of a symmet- 
ric triple well potential which diffuses symmetrically to 
the left and the right well. At a finite temperature and in 
absence of any bias force the particles are activated only 
by inherent thermal fluctuation resulting in equalization 
of population in the two side wells. However, if, in addi- 
tion, we allow the potential wells to roll by an external 
periodic signal, the escape over the potential barrier is 
modified by the interplay of the thermal fluctuations and 
coherent external signal, due to Stochastic Resonance^- 
9 . On the other hand when the height of the potential 
barrier is made to oscillate symmetrically or fluctuates 
around a mean value by the action of an external input 
signal the mean escape time over the fluctuating barrier 
exhibits a minimum at a particular value of frequency or 
correlation time of the external source due to Resonant 
Activation [TDHTB] . In both of these cases the time aver- 
aged escape rates from the middle well are equal and the 
stationary population of the left and right wells remain 
the same. 

Our aim of this letter is to explore a possible route 
leading to a splitting of the time averaged Kramers es- 
cape rate from the middle well due to the interference of 
these two resonances and to propose a convenient method 
for controlling the pathways of a parallel reaction for 
which the barrier heights corresponding to two product 
states are equal. For example, take the case of nucle- 
ophilic attack by X" (a halide ion of HX) at the carboxyl 
group of a ketone, say, Ri(R 2 )C = O producing D— 
R 1 {R 2 )C{OH)X and L— R 1 (R 2 )C(OH)X, two optical 
isomers (enantiomers) having same energy and stability 
but differing in their optical properties and hence bio- 
chemical activities. The middle well of the potential sig- 
nifies the reactant state and the terminal wells represent 
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the two product states of the parallel reaction. Specifi- 
cally, our objective here is twofold: First, to understand 
how the asymmetry in the time-averaged dynamics of the 
triple well potential driven simultaneously by two sinu- 
soidal forces results in differential average escape rates to 
two product states and unequal distribution of stationary 
population densities between them. Second, to explore 
the role of phase difference of the two interfering forces 
in determining asymmetric diffusion of the particles from 
the middle well and resulting localization in one state. In 
other words we look for a strategy for coherent control 
of pathways of a parallel reaction. As an interesting off- 
shoot of the analysis a selective process of enrichment 
of one of the two isoenergetic isomers under appropriate 
thermal condition is also explored. 

To illustrate the basic idea we begin with an over- 
damped Brownian particle in a symmetric triple- well po- 
tential V(x) (Fig.l) kept in a thermal bath at temper- 
ature T and subjected to two sinusoidal signals a\(t) = 
A\ sin (uiit + 4>\) and 02(f) = A% sin (u>2t + fa). The gov- 
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FIG. 1: (Color online) A schematic illustration the two config- 
urations of the potential under simultaneous action of the two 
signals ai(t) and a2(t) for A<j> = cj>i ~ 4>2 — and synchronized 
frequencies uii = u>2- 
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erning Langevin equation is given by 

jx = —V'(x) + A\ sin {u\t + (f>i) + A 2 a;sin(w 2 i + <fe) 

where V(x) — x 2 (bx 2 — c) 2 ; b and c are the parameters 
of the potential and 7 is the dissipation constant. Ui and 
0i (i = 1,2) are the frequency and phase of the signals. 
Thermal fluctuation of the bath is modeled by zero mean 
((r(i)) = 0) and delta correlation of noise, (T(t)T(t')) = 
2DS(t — if), D being the strength of the thermal fluc- 
tuation and is given by D = kT/~/. Here the additive 
signal ai(i) rocks the potential wells sidewise, whereas 
the multiplicative signal a,2(t) sets a symmetric oscilla- 
tion of the barrier height around AVq(= 4c 3 /27b) with 
an amplitude ±A2X%/2 at ±x m , respectively ( since the 
fluctuation is space dependent the amplitude of fluctua- 
tion of the barrier height around AVb is ±A2(x 2 n — x\)f2 
for terminal wells to middle well), where ±Xb and ±x m 
are the coordinates of two barrier tops and two termi- 
nal potential minima, respectively. The two configura- 
tions of the potential under simultaneous action of the 
two signals a\(t) and a2(i) are schematically illustrated 
in Fig.l for A(f> = 4>i ~ 4>2 — and synchronized fre- 
quencies ll>i — ll>2- As shown the barrier height for the 
transition from middle to right well fluctuates with an 
amplitude ± (A2XI/2 + A\Xb) whereas for the middle to 
left well the amplitude of fluctuation of the barrier height 
is ± (A2XI/2 — A\xb)- If the external modulations ai(t) 
and a,2{t) are small and very slow implying AVb 3> A^,A\ 
and the Kramers escape time (l/V^) for the unperturbed 
system is much smaller than the time period of the exter- 
nal input signals (1/Vfc <C 2tt juj\ = 2tt/uj2) one may con- 
sider the expressions for the time-dependent transition 
rates from the middle to left and right wells as follows 

W/L { + \ _ „ „ VT ^ [ A 1 x b sin(u 1 t+<f> 1 ) - g sin (c^ 2 t+<fe) , 
vv M\ l ) — 'fe ex PL D J' 

W*(t) = r k exp [+ AlXb sin +^ si " ( " 2t+02) ] , 
respectively. Here r k = ex P [ — T5^] ^ s the transition 
rate from the middle well for the unperturbed system; 
loq, u>b are the frequencies corresponding to the potential 
minimum (xo) and barrier top (xb), respectively. Fol- 
lowing McNamara and Wiesenfeld [8] if we expand the 
exponential term of the time dependent transition rate 
and keep the leading terms up to second order, the time 
averaged transition rates for tjj\ — U2 and (f>i = 4>2 = 
are given by 

(W&(t))t = r*jl + ^(Aix b -^) 2 | (2) 

It is apparent from the above expressions that as a re- 
sult of interplay of two resonances the period averaged 
transition rates from the middle to left and right wells 
significantly differ from each other at very low tempera- 
ture and tend to equalize in the high temperature limit. 



In Fig. 2 we present two representative plots for the ratio 
of transition rate to right and left wells as a function of 
Ft^JrrpWature and compare the result with simulations in 
the adiabatic regime l/r^ <C 27r/cji = 2ir/u}2- 

We estimate the ratio W r (= (W$) /{W^}) of the tran- 
sition rate as a function of external driving frequencies 
(synchronized frequencies, U)\ = 0J2 ) by standard numer- 
ical simulation of the Langevin equation using Huen's 
algorithm. We allow 10,000 test particles to move from 
middle well in either direction and count the number 
of particles which arrive in the left well(ni) and right 
well (tir), to calculate the ratio of the transition as given 
by W r = (Wm) I '(Wm) — n R /n L . It is apparent from 
Fig. 2 that the numerical analysis matches fairly well with 
our theoretical result. As a result of fluctuation of the 
barrier heights for both left and right wells, the transition 
rates (W^) and (W^) exhibit resonant activation inde- 
pendently as expected (when the Kramers escape time 
coincides with 2ir/oj 2 )• Moreover as the amplitude of 
fluctuation of the barrier height corresponding to the 
right well is larger (see Fig.l), the ratio W r differs sig- 
nificantly from unity and exhibits a resonance when the 
latter is plotted as a function of the frequency under a 
phase matched condition 4>\ = 4>2- 

In Fig. 3(a) we plot this ratio as a function of frequency 
between the synchronized input signals for several values 
of temperature to exhibit this asymmetry in resonant ac- 
tivation due to differential transition rate. In Fig. 3(b) 
we plot the ratio of the transition rate as a function of 
the phase difference of the two input signals for different 
values of strength of the input signal(A2). It is observed 
that the transition rate towards the side wells are equal 
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FIG. 2: (Color online) A comparison between analytical re- 
sult( based on analytical expression (2,3)) and numerical sim- 
ulation plotting W r as a function of temperature for the pa- 
rameter set: Ac/) = 0, wi = uj 2 = 0.0005, A\ = 0.1, A 2 = 
0.2, b = 0.1, and c = 1.0 
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(W r = 1) for the phase difference <f>\ — (j>2 — tt/2. The 
ratio of the rates W r can be inverted by reversing the 
phase difference from 0i — cf>2 = to <f>x — 4>2 = t 1 "- There- 
fore by controlling the phase of the input signals a\{t) 
and a2 (i), it is possible to manipulate the transition of 
particles from the middle well to the product states and 
hence the course of the parallel reaction. 

The above analysis is based on kinetic considerations. 
It is also worthwhile to turn our attention to the time av- 
eraged distribution of the particles in the two wells and 
the associated aspects of localization 14] . This is in the 
spirit of stationary product distribution of a parallel re- 
action. In this case we allow the trajectories to evolve 
dynamically for long time starting from an arbitrary ini- 
tial position in the well. Numerical simulation results 
show that the residential time distribution (Fig. 4) bears 
a marked asymmetry corresponding to a stochastic local- 
ization of the 

particles in right well for (f>i — 02 — and in the left 
well for 4>i — (p2 = 7T- For <pi — (j>2 — tt/2, the distribution 
of x(t) over time is more or less even for both the wells. 
In view of the input signal synchronization, a qualitative 
interpretation of this type localization may be given as 
follows: So long as the force 0.2 (t) causing symmetric os- 
cillation of the barrier height attains its lower value, the 
tilting force cti(i) points to the right well, so that the par- 
ticle in the middle well move towards the right well very 
quickly. On the other hand as the tilting force points to 
the left d2(t) sets the barrier height at a larger value and 
consequently the particle in the middle well takes rela- 
tively larger time to speed up from middle to left well 
for the simultaneous action of the synchronized signals. 
The particles in the middle well therefore have a greater 
chance to cross the right-hand barrier. In the right well, 
the amplitude of oscillation of the barrier height is much 
larger and so the particle coming into the right well es- 
cape from it more quickly and on returning back to the 
middle well it has again two options to cross the barrier 
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FIG. 3: (Color online) (a) Variation of W T as function of 
frequency for different values of temperature T — 0.2(dashed 
line), T = 0.25 (Solid line) and T = 0.4 (dotted line) and for 
the parameter set: A(f> = 0,Ai = 0.1, A% = 0.2, & = 0.1, c = 
1.0. (b) W r vs phase difference A(j> plot for different values 
of A2 and for the same parameter set as Fig. 3(a) but for 
lji = CJ2 = 0.05 and T = 0.2. 



as it was initially. Thus the particle dynamically spends 
most of the time in between the middle and the right 
well. To proceed further we require a quantifier which 
measures the asymmetry in localization in the two wells. 
To this end we choose the mean position of the particle 
as a measure of asymmetry. For a symmetric distribu- 
tion mean position (x) = and for the localization of the 
particles to left or right well, value of mean position is 
negative or positive, respectively. In Fig. 5(a) we present 
the variation of mean position as a function of synchro- 
nized input signal 

frequencies. With increase of the input signal fre- 
quency the mean position gradually shift to a maximum 
positive value followed by a decrease to zero at high fre- 
quency In Fig. 5(b) we show the mean position as a func- 
tion of phase difference between the two input signals. 
The mean position is zero for Acf> — tt/2 and it departs 
from zero as the phase difference differs from 7r/2. For a 
phase matched condition, <pi = 2 the particle are local- 
ized in the right well while for a phase reversal A(f> = tt, 
localization takes place in the left well. 

Can temperature influence the product distribution of 
a parallel reaction at the steady state? This question is 
intimately related to the manipulation of incoherent con- 
dition rather than coherence in selecting and controlling 
the reaction pathways. To have a closer look into this 
aspect we examine the variation of mean position (x) 
with temperature with the help of a discrete three-state 
model for the triple well potential. Three states are de- 
noted by xq, ±x m for the symmetric unperturbed system 
corresponding to three minima. The diffusional motion 
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FIG. 4: (Color online) Residential time distribution for sev- 
eral values of phase difference (A</> = (upper panel), A(f> = 
tt/2 (middle panel), A0 = tt (lower panel)) between two in- 
put signals, for the parameter set: Ui = 0J2 — 0.0013, A\ = 
0.1, A 2 = 0.25, T = 0.155, b = 0.1 and c = 1.0 
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causes transitions between them and it is schematically 
presented as 



L 



M 



R 



k R , (W M ), (W M ) denote the time averaged rate 
of transition from left to middle well, right to middle 
well, middle to right and middle to left well, respectively. 
The number of particles in the three states at time t 
are denoted by til, n R and ul- The governing master 
equations for n, (i — L, R, M) read as 



dn L 
~dT 
dn R 

dt 
dn M 

dt 



-k L n L + (W M ) n M 
-k R n R + (W M ) n M 
k L n L + k R n R - {(Wm) 



(4) 
(5) 

(W M ))n M (6) 



At the steady state (n'x, = ti'm = n R = 0) the proba- 
bility of finding the particles at the three wells Pi (i = 
L, R, M) are P L = (W M )k R /P, P R = (W M )k L /P and 



Ml 

Pm = k L k R /P where P = k R (W M ) + k R k L -f k l \v M / 
The expression for the mean position is then given by 



k L (W^ 



(x) with temperature is therefore expected to show a 
bell-shaped curve. The departure of (x) from zero to- 
wards positive direction indicates the preferential distri- 
bution of the product in the right well. The numerical 
simulation of the variation of mean position (x) as a func- 
tion of temperature for synchronized input signals under 
phase matched condition as shown in the Fig. 6 corrobo- 
rates this assertion. 

In summary, we have shown that depending on their 
phase difference, an application of two synchronized sig- 
nals on a particle in a triple well potential may lead to a 
splitting of the time averaged Kramers escape rate due to 
an interference of stochastic resonance and resonant acti- 
vation. This allows as to realize a strategy for achieving 
a preferential product distribution in the steady state 
of a parallel reaction. The present analysis thus reveals 
that in stochastic energetics[15. can be utilized to con- 
trol kinetically the pathways of a chemical reaction by 
appropriate manipulation of coherence and/or inherent 
thermal condition. 



xP(x) dx = x m P R + x Q P M - x m P L 



(W M )k L - (W M )k R 



2c J k R (W M ) + k R k L + k L (W M ) 



(7) 



The above expression clearly shows the dependence of 
mean position and probability on four time-averaged rate 
constants. Furthermore if we assume that k R and 
do not differ significantly then(a;), in general, turns out 
to be positive since by Eq([2p| (W M ) is greater than 
(W M ). Keeping in view of the Arrhenius temperature 
dependence of the individual rate constants, the variation 
of 
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FIG. 5: (Color online) position((a;}) vs frequency plot for 
several values of temperature and for the parameter set: 
A(j> = 0, Ax = 0.1, A 2 = 0.25, b = 0.1, and c = 1.0. (b) 
Mean position((:r}) vs phase difference (A(j>) plot for different 
input signal frequencies and for the same parameter set as 
Fig.5(a) but for T = 0.155. 
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FIG. 6: (Color online) Mean position jx^, vs temperature plot 
for different values of amplitude of the synchronized external 
signals and for the parameter set A<f) —0, uj\ = a;i=0.0005, 
b=0.1, and c=0.6 (all the quantities are dimensionless) 
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